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Abstract We show how liquid migrates in sheared un- 
saturated granular media using a grain scale model for 
capillary bridges. Liquid is redistributed to neighboring 
contacts after rupture of individual capillary bridges 
leading to redistribution of liquid on large scales. The 
liquid profile evolution coincides with a recently de- 
veloped continuum description for liquid migration in 
shear bands. The velocity profiles which are linked to 
the migration of liquid as well as the density profiles of 
wet and dry granular media are studied. 

Keywords Wet granular matter • Contact dynamics 
simulations • Liquid bridge • Cohesion ■ Liquid 
migration 



1 Introduction 

How does liquid spread in wet granular media under 
shear? We know from everyday experience that mixing 
large amounts of powders and liquid is a challenging 
task. Liquid bridges appearing in wet granular matter 
exert cohesive forces onto the grains resulting in sub- 
stantial changes of the mechanical properties |L highly 
enhancing the difficulty in obtaining a homogenous mix- 
ture of liquid and powder. The development of advanced 
techniques for the improvement of mixing processes re- 
quires gaining knowledge about the influence of the 
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grain scale redistribution of liquid onto the large scale 
liquid migration in a granular medium. 

The fundamental difference between ordinary cohe- 
sive and wet granular matter is that in the wet case 
the cohesive force acts only after liquid bridge forma- 
tion exhibiting hysteresis, resulting in a rich variety of 
observed phenomena [TJ[2]. The interaction distance of 
two grains is determined by the rupture distance of 
a capillary bridge. The forces and rupture distances 
were experimentally measured in Ref. 3, where the au- 
thors provided well established formulas suitable for 
numerical simulations. Discrete clement methods are 
widely used to model wet granular matter. The sim- 
plest model to account for hysteresis is the minimal 
capillary model [1 . Simplified models for capillary clus- 
ters were introduced in Ref. :4 . Most simulations imple- 
ment geometry imposed or homogenous capillary bridge 
volumes [5]. Richefeu [5] implemented a redistribution 
scheme for liquid and investigated the mechanical prop- 
erties of wet granular matter in direct shear but did 
not allow for larger changes in liquid concentrations. 
A very important problem which was rarely studied in 
the past is the transport of liquid in wet granular mat- 
ter. It is still unclear which mechanisms are responsible 
for an altering liquid concentration. M. Scheel et al. [5] 
studied the equilibration of liquid bridges through thin 
films on the particles due to Laplace pressure differ- 
ences. However, liquid is also redistributed due to rel- 
ative motion of grains leading to steady rupture and 
creation of liquid bridges. In this paper, we study this 
redistribution process in detail and present a model for 
capillary bridges taking into account the redistribution 
of liquid after bridge rupture. As opposed to previous 
studies where the volume of the bridges is either ho- 
mogeneous or given by the local geometry we explicitly 
account for the volume of individual capillary bridges 
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given by the amount of liquid which was trapped by 
the roughness of the particles. The redistribution of liq- 
uid gives rise to locally changing liquid concentrations 
which will be studied in periodic simple shear between 
two rough walls [7J|8] . 

2 Model 

We model wet granular matter at low liquid content 
under shear. For this purpose we use Contact Dynamics 
to model spherical hard particles and develop a model 
for capillary bridges with adhesive forces. 

2.1 Particle dynamics 

Particles are modeled using Contact Dynamics origi- 
nally developed by J. J. Moreau which is suitable for 
rigid particles [51IT0]. Here, contact forces are calcu- 
lated based on perfect volume exclusion and Coulom- 
bian friction. In the framework of Contact Dynamics, 
more contact laws have been established such as rolling 
friction, hydrodynamic lubrication and cohesion |11[|121 
[T5] . The fundamental difference to standard discrete el- 
ement methods lies in the determination of the contact 
forces which are calculated based on perfect volume ex- 
clusion of the contacting particles. For instance, the 
contact normal force F n between two grains is calcu- 
lated such that the two grains do not overlap at the 
next time step. A tangential constraint force F^ est is 
subsequently calculated under the condition of no slip 
at the contact which requires the relative tangential 
velocity to vanish at the next timestep. However, |F t | 
must not exceed the threshold value fi\F n \ where fi is 
the friction coefficient. If the previous calculation of 
latest | - g g rea ter than the threshold we set \F t \ = fJ,\F n \ 
which is a sliding, energy dissipating contact, else we 
set F t = F t test . 



the liquid air interface and s the separation [3 . This 
empirical formula for the forces is well established and 
valid for small capillary bridges up to V « 0.03r 3 . 
The contact angle is set to zero in this work. We 
see that the capillary force is independent of the vol- 
ume when s — 0. However, roughness influences the 
attractive force as soon as the length scale of the liquid 
bridge is comparable to the length scale of the rough- 
ness [T3]. In that case, the force can be much smaller. 
The authors distinguish three regimes for very small, 
for intermediate and for large amounts of liquid. The 
first two regimes are characterized by a monotonically 
increasing capillary force with volume before the force 
becomes constant in the third regime. We simplify the 
force law by choosing a threshold value below which the 
attraction vanishes. 




Fig. 1 (Color online) The liquid redistribution scheme after 
bridge rupture. The total amount of liquid is conserved in 
this model. 



2.2 Capillary bridges 

A wetting liquid forms a wetting layer on a grain sur- 
face. As soon as two grains touch, liquid from the wet- 
ting layer accumulates at the contact point to form a 
capillary bridge. This meniscus gives rise to an adhe- 
sive capillary force. The force was measured in detail in 
Ref. [3] and the authors derived an empirical formula 
for the capillary force given by 



1 + l.Obsy^JV +2.5s 2 r/V 

where V is the bridge volume, r the average radius of 
curvature of the two spheres, r the surface tension of 



When two grains separate, the liquid bridge eventu- 
ally ruptures. The above mentioned authors also pro- 
vide an empirical formula for the rupture distance s c 
which is 

s c = V YjZ (2) 

for zero contact angle and small volume. 

There are two relevant mechanisms leading to a re- 
distribution of liquid. First, in case of rupture, the liq- 
uid is sucked back onto the grains very fast. At low 
water contents, the pressure in the bridges compared to 
the pressure in the wetting layers is small, such that the 
liquid is drawn into the already existing bridges 1 . The 
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time scale for this redistribution process is expected to 
be of the same order as the time scale for bridge for- 
mation between two grains when the roughness is satu- 
rated with liquid. If the relative velocities of the grains 
are small enough this redistribution of liquid can be 
considered to take place instantaneously. 

Secondly, as shown experimentally in Ref. [15] , there 
is a flux of liquid between capillary bridges through the 
wetting layers on the bead surfaces or through the vapor 
phase driven by Laplace pressure differences. However, 
the equilibration of capillary bridges can be neglected if 
the lifetime of a capillary bridge is much smaller than 
the equilibration time scale. For a fixed equilibration 
time scale, this criterion is fulfilled for sufficiently high 
shear rates since the relative velocities of the grains 
become larger and the lifetimes of the capillary bridges 
become smaller with increasing shear rate. 

Therefore, a model in which the equilibration of liq- 
uid bridges is neglected and in which the liquid is in- 
stantaneously redistributed to neighboring bridges after 
bridge rupture is indeed valid in a certain range of shear 
rates: The relative motion of the grains must be slower 
than the redistribution of liquid after bridge rupture 
but faster than the equilibration of liquid bridges. The 
time scale for bridge formation is less than a millisec- 
ond for grains of one millimeter diameter [TinBj . On the 
other hand, typical equilibration times of liquid bridges 
are of the order one to five minutes [TS] for glass beads of 
half a millimeter diameter which is a much larger time 
scale. Since we expect the time scale for liquid redistri- 
bution after bridge rupture to be of the same order as 
the time scale for bridge formation, the assumption of 
instantaneous liquid redistribution after bridge rupture 
and neglection of equilibration of liquid bridges will be 
valid in the range 10 _2 s _1 <C 7 <C 10 3 s~ 1 which spans 
a large range of applicable shear rates 7. 

Based on these considerations we propose the fol- 
lowing model for redistribution: Particles can carry an 
arbitrary amount of liquid Vf > V m in in a wetting layer 
where V m i n is a fixed amount of liquid which is trapped 
in the roughness [2] . When two particles come into con- 
tact and at least one particle has Vf > V m in , a capillary 
bridge is instantaneously formed. Since for each grain, 
V m in is fixed, the available liquid for bridge formation 
is Vf — V m in therefore, without loss of generality, we can 
use Vmin — and the bridge volume V is given by the 
sum of the two involved films, i.e. after bridge formation 
V = Vj + Vf and Vf = 0, k e 1,2. If the separation 
of particles exceeds the critical distance s c , the bridge 
volume V rup is equally split between the two particles 
Pj such that Vj — V rup /2, j £ {1,2}: At the same 
time, if there exist further contacts, namely capillary 
bridges or dry contacts, all liquid is equally distributed 



among them such that the new volume of a neighbor- 
ing contact Cf, i G {1 ■ ■ • Nj} of particle Pj is given by 
V?> new = V rup /(2N j ) + Vf where Nj is the num- 
ber of neighboring contacts of particle j. An example 
of this redistribution scheme is shown in fig. [1] where 
initially, all bridges have the same volume. If Nj = 
the liquid is kept in the wetting layer of particle j for 
the formation of a bridge at a later time. To prevent 
formation of liquid clusters via bridge coalescence [5] 
we use an upper threshold V ma x of the bridge volumes 
which must not be exceeded. In that case, all contacts 
of that particle with V < V m ax are filled first and if 
all are filled completely, the liquid remains in the film. 
This case however, happens very rarely since we con- 
sider relatively low liquid contents only. In our model, 
liquid mass conservation is ensured, whereas in suction 
controlled models, e.g. [H| the Laplace pressure is pre- 
scribed such that the total amount of liquid can change. 



3 Simulated Setup 

We simulate normal stress controlled periodic simple 
shear between two rough walls. We use particle radii 
uniformly distributed between 0.8 and 1 in units of 
the largest particle radius R. Fig. [2] shows the geom- 
etry used in our numerical experiments. The bottom 
wall is fixed at position z = but moving at con- 
stant speed v s hear along the x direction. The top wall 
is fixed in x direction and a pressure P is applied onto 
it such that the position of the wall z wa u(t) fluctuates 
around z wa ii — (z wa ii{t)). Gravity is neglected and we 
use periodic boundary conditions in x and y direction. 
The system dimensions are L x = 2QR, L y = 12R, 
Zwaii ~ 78R. The results will be presented in terms 
of the global shear strain 7, which is the shear dis- 
placement v s heart divided by the system height L z . 
As proposed by Rognon et al. [T7] there are two di- 
mensionless parameters which govern the mechanics of 
the system, the inertial number I = jR^/ p/P where 
7 = dv x /dz is the local shear rate and the cohesion 
number 77 = F r c nax /PR 2 = 2-KPjPR which is the ra- 
tio of the largest cohesive force to the force exerted on 
a particle by the pressure P. We first study the sim- 
plest case 77 = where shear flow is approximately 
homogeneous which can be realized in experiments by 
choosing a sufficiently high pressure. We subsequently 
study the case where r\ > and in both cases, we 
fix / rs 0.008. Initially, we impose a Gaussian distri- 
bution of liquid bridge volumes. As a function of the 
bridge position z the volumes are initialized according 
to Vb(z) = Aexp(— (z — z wa u/2) 2 /(Jq) with amplitude 
A and width ctq and we monitor the evolution of the liq- 
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uid distribution during shear. Here, the shearing walls 
are assumed to be hydrophobic. 

4 Results: r) — 

In this section, we present results in the limit of van- 
ishing cohesion number. We show the obtained velocity 
profiles and the consequences for the liquid distribution. 

4.1 Velocity profile 

Velocity profiles have already been studied for non- 
cohesive systems in detail [71IT5] . It was shown that un- 
der high confining pressures the system is alternating 
between diffuse shear and shear banding occurring at 
random height z. These shear bands are however not 
stable, in fact they disappear and reappear constantly 
at different positions. However, when averaging the ve- 
locity profile over long time scales it becomes linear 
apart from a slight S-shape which would vanish in the 
limit / — > [TB] . The time averaged velocity profile for 
our system is shown in the inset on the left of fig. [3] 

4.2 Liquid profile 

During shear the liquid profile changes. In fig. [2] we ob- 
serve that the initial Gaussian liquid distribution (left 
picture) spreads towards the top and bottom walls (mid- 
dle and right pictures) . There are two relevant processes 
involved causing the spreading of liquid. It is known 
that in plane shear flows particles undergo a diffusive 
motion and therefore, also liquid which is carried by 
the menisci will diffuse in space QI5]. Secondly, there is 
a transport of liquid associated to liquid bridge rupture. 
As explained in section |2"T21 the liquid is redistributed to 
all neighboring liquid bridges after bridge rupture which 
means that locally, after a bridge rupture event, there is 
a liquid flux away from the rupture point. Fig. [3] shows 
the evolution of the liquid distribution at shear strains 
7 = 0, 30, 60 where Q is the average bridge volume per 
particle. The data was averaged over a small strain in- 
terval, five independent runs and in slices of width 2R. 
We observe the liquid distribution to remain Gaussian 
(solid lines) with an increasing width (7(7) which is ex- 
actly the case for diffusive processes. The solution of 
the standard diffusion equation with Gaussian initial 
condition is a Gaussian for which a 2 (t) increases lin- 
early in time |20j . Since 7 is proportional to t it can be 
expressed as a 2 (7) = a 2 + 4Dj where D is the diffusion 
coefficient. The inset in fig. [3] shows a 2 as a function of 
7 (blue crosses) and an excellent linear fit to the data 




Fig. 2 (Color online) Snapshots of plane shear between two 
walls at the beginning (left), at an intermediate time (middle) 
and at the end of the simulation (right). Only the central 
part of the sample is shown and the y direction points into 
the plane. The capillary bridges are indicated by dark blue 
lines connecting the centers of two spheres whose width is 
proportional to the bridge volume. It is seen that the liquid 
is spreading towards the top and the bottom wall. 



(red line) indicating diffusive behavior. Here, D is about 
0.4i? 2 per unit shear strain. Recently, a continuum de- 
scription was developed for a similar model which was 
experimentally verified in a split bottom shear cell [21] ■ 
After bridge rupture, liquid was redistributed to neigh- 
boring capillary bridges taking into account the Laplace 
pressure in the bridges as well as the distance from the 
rupture point to the neighboring bridges. A modified 
diffusion equation was derived by considering that liq- 
uid fluxes happen via bridge rupture events which are 
proportional to the volumes of the ruptured bridges and 
to the local rupture rate. The continuum equation was 
given by 

Q b = C-^(BQ b ) (3) 

where Qb is the average bridge volume, B is the bridge 
rupture rate and C is a constant. We notice that the 
equation is applicable to the plane shear geometry as 
well: Since the velocity profile is linear, we expect the 
bridge rupture rate to be constant such that we recover 
the ordinary diffusion equation Qb = Dd 2 /dz 2 Qi>. As- 
suming that the amount of contacts per particle is con- 
stant, Qb can be replaced by Q. 



4.3 Dependence on liquid content 

The influence of the liquid content on the diffusion con- 
stants was investigated by varying the amplitude A 
of the Gaussian liquid distribution. Here, we used a 
smaller system for computational reasons. Fig.2]shows 
the diffusion coefficients D for different amplitudes A 
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Fig. 3 (Color online) The averaged liquid content as a func- 
tion of the height z in periodic simple shear for r\ = at 
strains 7 = 0, 30, 60. The lines are Gaussian fits to the data 
points. The inset on the right shows a linear increase of the 
width a 2 with time as expected for diffusive behavior. The 
inset on the left shows the averaged velocity profile as a func- 
tion of the height z. 




0.02 ^4 0.03 



0.05 



Fig. 4 Diffusion constants D as function of the amplitude 
A of the initial Gaussian liquid distribution. They decrease 
with increasing amplitude because the rupture distances are 
larger and thus the rupture rate decreases. The solid line is a 
fit according to D(A) = o/(l + bA 1 / 3 ). 



of the Gaussian. It can be seen that the diffusion coeffi- 
cients decrease with increasing liquid content. A certain 
dependence is indeed expected: Since the rupture dis- 
tance of a capillary bridge scales as s c ~ the strain 
needed to rupture a capillary bridge increases with in- 
creasing bridge volume. The strain needed to break a 
contact in a dry granular medium is of order unity |17) 
such that the rupture time for dry granular media can 
be estimated by T Sc= q ~ I/7. The additional strain j r 
needed to rupture a capillary bridge can be estimated 
by 7r ~ Sc/ R- Thus, the rupture time T Sc —q for dry con- 
tacts is increased by T Sc oc T Sa=0 + s c /Rj oc 1 + s c /R 
such that we expect the diffusion constants to decrease 
with Q as D ~ 1/(1 + bQ 1 / 3 ), where b is a fit parame- 
ter. Fig. [4] shows the diffusion coefficient D for different 
amplitudes of the Gaussian which can be well fitted by 
D(A) = a/ (I + bA 1 / 3 ) with fit parameters a and b. 



5 Results: r) > 

Now we study the case of finite cohesion number rj. We 
obtain a non-homogenous velocity profile and investi- 
gate its influence on the liquid migration. 



5.1 The non-homogenous case: velocity and density 
profiles 

As opposed to the case rj — here the capillary forces 
play a dominant role since weaker tangential stresses 
are less capable to rupture the bonds between the grains. 



Since wet grains stick together, the velocity profile de- 
velops a plateau and the shear rate drops inside the 
wet area. Examples for such plateaus can be seen in 
fig. O for two different cohesion numbers. The shearing 
of a mixture of two granular materials with different 
friction coefficients was studied recently by Unger [22] . 
The upper half of the shear cell was filled with material 
having a larger coefficient of friction than that of the 
lower part. He observed that in such a system, the up- 
per part with larger friction coefficient undergoes shear 
hardening. The shear rate was found to be given by 

7 oc cosh(k(z — z c )) (4) 

where z c is the middle position of the upper part of 
the cell and therefore, the velocity profile is given by a 
sinh function. Since in our case friction is effectively en- 
hanced due to the attractive force, the same arguments 
as in Ref. [22] are applicable to our system. Indeed, 
fig. [5] shows typical velocity profiles appearing during 
shear which can be fitted by E + Asmh(k(z — z c )/z wa ii) 
inside the wet region. Since the liquid profile changes, 
we only show velocity profiles averaged over a small 
strain interval of 0.4. We see that increasing the cohe- 
sion number 77 leads to a non linearity in the velocity 
profile which becomes more pronounced with increas- 
ing rj. The parameters E, A, k and z c fluctuate very 
much over time as shown exemplarily for the param- 
eter k in fig. [S] We observe k to significantly increase 
with increasing cohesion number 77. As pointed out by 
Unger, the amount of agitation a layer in the x-y plane 
receives from the next upper respectively lower layer 
decreases with increasing k. Here, the same behavior is 
found, i.e. the shear rate in the wet region drops more 
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Fig. 5 Typical velocity profiles emerging during shear. Av- 
erages were performed over five time steps during a strain 
of about 0.4. The plateau in the velocity profile is flatter for 
larger r\. Lines are fits to the data. 
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Fig. 6 (Color online) The coefficient k for different cohesion 
numbers r\ as a function of the strain 7. For increasing 77, 
the coefficient k increases which reflects the hardening of the 
granular medium inside the wet region. 



for larger cohesion numbers. The explanation for the 
appearance of a velocity profile given by a sinh func- 
tion in Ref. [22] was based on the presence of a shear 
stress heterogeneity. However, we could not find notice- 
able stress heterogeneities in our systems such that we 
believe that in our case, a larger effective viscosity of 
the granular medium caused by cohesion is responsible 
for the locally decreased shear rate. 

Next, we study the density profiles for different co- 
hesion numbers. The density 4> of cohesive granular me- 
dia in plane shear was found to decrease with increasing 
77 [17] , due to the formation of large stable pores in the 
medium during shear. Somewhat counterintuitively, we 
find that the density shown in fig. UJ for three different 
cohesion numbers 77 increases in the wet region despite 
the presence of cohesive forces. On the other hand, the 
shear rate respectively the local inertial number in the 
wet region is decreased which in fact inhibits the for- 



Fig. 7 (Color online) The density <f> as a function of the 
height z for different cohesion numbers rj. We see that the 
density increases inside the wet region when increasing the 
cohesion number rj. 



mation of larger pores. Our observation coincides with 
Ref. [22] where density differences originate from hav- 
ing two materials with different friction coefficients, al- 
though the microscopic origin for increased density is 
different in our work. 



5.2 Liquid profile: The non-homogenous case 

We showed that under shear the liquid migrates due to 
diffusion and liquid bridge rupture. Both processes are 
however, a function of the shear rate. Since the shear 
rate drops in the wet region, although the global shear 
rate is the same as for the case ry = 0, we expect the liq- 
uid to migrate on larger time scales which reflects a fre- 
quently encountered problem in industrial applications 
when mixing granular matter with liquid. If the humid- 
ity profile was non-homogenous but the grains were wet 
everywhere, we would recover the linear velocity profile 
again which enhances the mixing properties. We can al- 
ready deduce that mixing dry and wet powder is even 
more difficult than equalizing the liquid distribution of 
completely wet samples because of cohesion and the ap- 
pearance of powder clumps. The evolution of the liquid 
profile for rj = 0.47 is shown in fig. [BJ As opposed to the 
case 77 = we see a deviation from a Gaussian at the 
border between the wet and the dry region at z ~ 50. 
It corresponds to the place where the shear rate pro- 
file starts to rapidly deviate from a constant. At this 
point, we cannot describe anymore liquid migration as 
a purely diffusive picture but it is necessary to take the 
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Fig. 8 (Color online) The averaged liquid content as a func- 
tion of the height z in periodic simple shear for r) = at 
strains 7 = 0, 80. The dashed lines are Gaussian fits to the 
data points and the black solid lines are the numerical solu- 
tions of eq. (J5J. The inset on the left shows the liquid content 
in different slices as a function of the strain where the solid 
lines are results from simulations and the dashed lines are the 
solutions of eq. (|SJ. The inset on the right shows In Q as a 
function of z, where filled diamonds are results for r\ = 0.47, 
crosses for r\ = and dashed lines are Gaussians. 

shear rate profile into account. We apply the continuum 
description 

Q = C-^(iQ) (5) 

from Ref. [21] to our system. To solve eq. (O numeri- 
cally we use the measured shear rate profile as an input 
to eq. ([5]) and the constant C was set to 0.475 in units of 
R 2 . The solid line in fig. [8] shows the numerical solution 
of eq. ([5]) which corrects the deviations from a Gaus- 
sian at the interface between dry and wet granulate. To 
have a better view on the accuracy of eq. ([3]) we plot 
ln(Q) as a function of the height in the inset on the 
right (filled symbols) in fig. [5] as well as the results for 
r\ = for comparison. We see that the deviations from a 
Gaussian are much more pronounced for i] = 0.47 than 
for r\ — caused by the shear rate. This shows that in 
our model, the shear rate profile determines how liq- 
uid migrates within the system and must be taken into 
account in a diffusive description for liquid migration. 
The reason why we see deviations from a Gaussian for 
?y = is on the one hand due to the boundaries and on 
the other hand due to the slight non-constancy of the 
shear rate profile. The observation of purely diffusive 
liquid migration in a numerical experiment would be 
possible in larger systems and by choosing appropriate 
boundary conditions without walls, e.g. Lees-Edwards 
boundary conditions [23] . The inset in fig. [8] on the left 
shows the evolution of Q in slices corresponding to the 
data points in the main panel as a function of strain (red 
solid lines) which is followed by the numerical solution 



(black dashed lines) of eq. §5§ very nicely. Since the local 
shear rate drives the liquid migration, improving mix- 
ing properties of grains and liquids can be achieved by 
decreasing the cohesion number. Experimentally, this 
requires increasing the pressure onto the top wall. 

Since the liquid migrates during shear, the velocity 
profiles are expected to change in time. The shear rate 
can be fitted by j(z) = c+a(z — z c ) 2 with fit parameters 
a, c and z c for sufficiently small \z — z c \. In our simula- 
tions, we found that c does not change much with strain 
for fixed rj. Therefore, the parameter a is a measure for 
the width of the plateau which increases with decreas- 
ing a. Fig. [9] shows the parameter a as a function of 
7 for two different cohesion numbers. In both cases, a 
decreases consistently with strain. The fact that the ve- 
locity profile is affected by the liquid distribution could 
be used in experiments to indirectly measure the posi- 
tion of the interface between the dry and wet region by 
measuring the velocity profile of the granulate. 



x 10 




20 40 7 60 80 100 

Fig. 9 The parameter a of the fit of the shear rate is plotted 
as a function of the strain 7. Since liquid is migrating, the 
plateau in the velocity profile widens which is reflected by 
the decreasing fit parameter a with 7. 



6 Conclusion 

We studied liquid migration in unsaturated granular 
media in plane shear. We found that liquid migrates 
diffusively if the shear rate profile is constant along the 
sample in agreement with experimental observations 
made in process engineering [24U25] . However, a non- 
homogeneous shear rate profile leads to a different liq- 
uid migration pattern in agreement with Ref. |21) . The 
influence of the cohesion number 77 on velocity profiles 
was studied and linked to the evolution of the liquid 
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profiles. Homogenous mixtures of liquid and powder are 
achieved by increasing shear rates and stresses. Direct 
experimental measurements of the liquid distribution in 
sheared granulates are difficult, but would be feasible 
in indirect measurements of the velocity profiles since 
the liquid distribution affects the velocity profile. 
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